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Introduction 

In the past year this group continued essentially full time research on extraterrestrial 
materials, and the question of the origin of the solar system. The continuing scientific staff consists 
of the P.I. and Visiting Scientist Michael Petaev. Vitae for Wood and Petaev appear in Sec. 6. We 
benefit from the part time services of a Project Administrator (Judith Terry) and a Secretary 
(Muazzez Lohmiller). 

In January 1999 the P.I. assumed the Chairmanship of COMPLEX, the Committee on 
Planetary and Lunar Exploration of the Space Studies Board, National Research Council. 

Wood and Petaev were authors or coauthors of 21 publications, new manuscripts, and 
abstracts in the last year. These are listed above, and referenced by number [n] in the discussion 
below. Other references to the literature made in this Section are listed in Sec. 3. 


New PCMET and CWPI condensation codes 

Two new Fortran codes - PCMET (for Co-bearing systems) and CWPI (for Mn-bearing 
systems) - have been developed from our previous codes PHEQ (Wood and Hashimoto, 1993) and 
PHEQ4 (Petaev and Wood, 1998a) which model nebular condensation in the H-C-O-Na-Mg-Al-Si- 
S-Ca-Fe and H-C-N- 0 -Na-Mg-Al-Si-P-S-Cl-K-Ca-Ti-Cr-Mn(Co)-Fe-Ni systems, respectively. Both 
codes implement the CWPI condensation model (Petaev and Wood, 1998), which allows for 
programmed withdrawal of minerals that are no longer accessible to reaction as they condense. Solid 
solution models and the thermodynamic database were substantially modified to resolve 
inconsistencies reported in our previous work (Petaev and Wood, 1998b). The new solid solution 
models that have been added are as follows: corundum (A1 2 0 3 - Cr 2 0 3 ), olivine (Mg 2 Si0 4 - Fe 2 Si0 4 

- CaMgSi0 4 - CaFeSi0 4 - Ni 2 Si0 4 ), orthopyroxene (Mg 2 Si 2 0 6 - Fe 2 Si 2 0 6 ), melilite (Ca 2 Al 2 Si0 7 - 
Ca 2 MgSi 2 0 7 ), plagioclase (CaAl 2 Si 2 0 8 - NaAlSi 3 0 8 - KAlSi 3 O g ), clinopyroxene (CaMgSi 2 0 6 - 
CaAl 2 Si0 6 - CaTi 0 sMg 0 5 AlSi0 6 ), metal (Fe - Ni - Co - Si - Cr, Co ideal), cordierite (Mg 2 Al 4 Si 5 0i 8 

- Fe 2 Al 4 Si 5 0 18 ), and spinel (MgAl 2 0 4 - FeAl 2 0 4 - FeCr 2 0 4 - MgCr 2 0 4 ). Sources of thermodynamic 

data are listed in [1]. 

To test our codes we compare our condensation sequence calculated for a nebula of solar 
composition (Anders and Grevesse, 1989) at P tot =10' 3 bar with that of Yoneda and Grossman (1995) 
[YG95 ], in Table 1. It should be noted that our database of condensed phases is ~5 times larger than 
that of [YG95]. Our condensation sequence includes all the phases of [YG95] plus Ca 4 Ti 3 Oi 0 , 
Ca 3 Ti 2 0 7 , barringerite, and schreibersite, which were not included in their database. In order to 
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Table 1. Condensation temperatures (°K) of phases facilitate comparison between our condensation 

sequence and that of [YG95] we removed 
Ca 4 Ti 3 Oio, Ca 3 Ti 2 0 7 , and P-bearing compounds 
from our database and calculated the new 
condensation sequence listed in Table 1 . Our 
condensation temperatures for corundum, 
perovskite, melilite, Al-spinel, Cr-spinel, metal, 
olivine, orthopyroxene, Ti 3 0 ; -, and TL 1 O 7 are 
consistent with those of [YG95] within 1 1°K, 
which we consider to be excellent agreement 
considering that there may be minor differences 
in the thermodynamic data between [YG95] and 
us. The major discrepancy in Table 1 is in the 
condensation temperatures and compositions of 
clinopyroxene and plagioclase solid solutions. In 
our condensation sequence clinopyroxene 
condenses earlier by 23°K and is initially Ti, Al- 
rich, with rapidly increasing Di content; in 
contrast, the initial clinopyroxene of [YG95] is 
pure diopside. This discrepancy is undoubtedly 
due to the ideal solid solution model of clinopyroxene used by [YG95]. Moreover, our composition 
of Ti, Al-rich clinopyroxene is fully consistent with the compositional zoning in the clinopyroxene in 
Fluffy Type A CAIs (e.g., Brearley and Jones, 1998) which are considered to be nebular condensates. 
The greater stability of clinopyroxene in our database is unquestionably responsible for the large 
difference in condensation temperatures of sphene in Table 1. The difference in plagioclase 
condensation temperatures is probably caused by differences in the solid solution models. There are 
also substantial differences in the stability fields of hibonite and grossite between [YG95] and our 
work which result from differences in thermodynamic data used for these minerals. We prefer the 
more’precise data of Kumar and Kay (1985) which predict stability fields of hibonite and grossite 
that are consistent with the earlier calculations of Kornacki and Fegley (1984). 

The PCMET/CWPI code, run on a Windows NT-based PC (333 MHz processor, 64 MB 
RAM) typically needs 20-30 sec to calculate an equilibrium mineral assemblage. At the 31 s LPSC 
we made our source codes and databases available to any interested party. 

Thermodynamic model of condensation of Bencubbin/CH-like chondrite metal 

Our earlier work on modeling condensation of CH chondrite metal [2-5] has been extended to 
Bencubbin/CH-like chondrites (QUE9441 1 and HH237), which are petrographically and chemically 
similar to CH chondrites but have higher metal abundances and contain more (15-20 vol. %) large 
zoned metal grains (Weisberg et al„ 1999; Krot et al., 2000). Compositional profiles across eight 
large zoned metal grains of the QUE9441 1 were measured by our collaborators A. Meibom and A. 

N. Krot at the University of Hawaii. The profiles (see appended manuscript) show Ni concentrations 
gradually decreasing from cores of the grains to their rims. Co generally follows Ni, but many 
profiles are somewhat ragged. Most Cr profiles show fairly uniform Cr concentrations in the cores 
with gradual increase in the peripheral portions of the grains. In order to reproduce the 
compositional profiles in these grains we carried out detailed modeling of non-equilibrium metal 
condensation in nebular systems with non-solar dust/gas ratios [6,7] using our new PCMET code [1]. 
Calculations were made at a nebular pressure of 10^ bar and for several sets of three parameters - Cr 
depletion factor, isolation degree q, and system dust/gas ratio - which were estimated from the 


stable in a system of solar composition at P t ot=10 bar. 
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mineralogy and chemistry of the QUE9441 1 and HH237 chondrites. Details of our model are 
described in the appended manuscript [7]. 

It was found that the condensation of these grains took place under conditions that differed 
for each grain, indicative of a high degree of heterogeneity in the formation region of Bencubbin/CH- 
like chondrites. The relatively high £, values (0.3-5 rel.% per K of cooling) required for condensation 
of these grains are attributable to non-equilibrium condensation resulting from rapid growth of the 
zoned metal grains. It was shown that the raggedness of the Co zoning profiles is apparently due to 
the growth of grain cores by coagulation of small, compositionally diverse grains of early Ni,Co-nch 
condensate rather than by direct deposition of Fe, Ni, Co, and Cr from the cooling nebular gas. The 
smooth compositional zoning of all elements in the grain rims is consistent with their direct 
deposition of Fe, Ni, Co, and Cr from the cooling nebular gas onto grain cores. Finally, zoned metal 
grains of the Bencubbin/CH-like chondrites must have been transported from their formation regions 
in order to escape secondary low-temperature alteration. 


CR clan chondritic materials: Further support for the CWPI-type condensation 

The recently recognized CR clan of meteorites includes the CR, CH, and Bencubbin/CH-like 
chondrites and the Acfer 182 and LEW 85332 ungrouped chondrites (Weisberg et al., 1995; 

Weisberg et al. 1999, Krot et al., 2000). These meteorites are quite different petrographically, but 
their similar mineralogical, chemical, and isotopic characteristics imply that they must have formed 
in the same nebular region from closely related reservoirs (Weisberg et al., 1995). Except for the 
original CR chondrites, other members of the CR chondritic clan have largely escaped secondary 
alteration, so they provide valuable information on the physicochemical conditions in their nebular 
source regions. During the last year we continued collaborative studies of these meteorites. 

The CH chondrites contain a population of small (15-30 pm) ferrous silicate spherules which 
consist of cryptocrystalline ol-opx-normative material, SiCVrich glass, and rounded-to-euhedral 
FeNi metal grains. The silicate portions of the spherules are highly depleted in refractory lithophiles 
(Ca, Al,Ti) and enriched in SiOi, FeO, MnO, Cr 203 , and Na20 relative to the dominant magnesian 
chondrules of the CH chondrites. The Fe/(Fe+Mg) ratios in silicates show strong positive correlation 
with Fe concentrations in coexisting metal grains, implying that this correlation is not due to a redox 
exchange of Fe between silicate melt and metal but is instead a nebular feature. We used our PCMET 
[1] and METEOMOD (Ariskin et al., 1997) thermodynamic codes to model the nebular condensation 
of solid precursor materials of these spherules, and their subsequent melting in a chondrule-forming 
event. Metal liquidus temperatures were calculated from an assessment of experimental data on the 
Fe-Ni system (Lee, 1993). It was found [10,1 1] that the observed properties of these spherules can be 
explained by non-equilibrium condensation of solid metal and silicate precursors in dust-enriched 
nebular systems (up to 300x solar) at 10 4 bar, followed by remelting of the precursor at relatively 
high f 02 (2. 1-3.1 log units below IW buffer) and crystallization of metal from supercooled (13- 
1 1 1°K) silicate melts. 

Many Mg-rich porphyritic chondrules (Type I) of the CR chondrites are surrounded by Si0 2 - 
rich igneous rims which consist of chemically zoned Mn- and Cr-rich pyroxenes, glassy mesostases, 
rounded FeNi metal blebs, and a silica phase. The mesostases are enriched in Na and K and depleted 
in Ca relative to those of the host chondrules. [12] concluded that the chemical fractionations among 
precursor materials of the host chondrules and their rims were produced by isolation and subsequent 
melting of earlier, high temperature condensates in the course of nebular condensation in the 
formation region of CR chondrites, in accordance with condensation sequences predicted by the 
CWPI condensation model of Petaev and Wood (1998a). 
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Secondary alteration in carbonaceous chondrites 

Petaev continued collaboration with the Hawaii group in studies of aqueous alteration in 
carbonaceous chondrites, which resulted in series of publications [13-16]. He carried out detailed 
thermodynamic analysis of the stability of the metal-sulfide-magnetite-Fe,Ni-carbide-fayalite 
assemblages observed by his collaborators in the reduced CV chondrite breccia Vigarano. It was 
found that carbidization of the metal might have taken place in either a nebular or an asteroidal 
environment (in the absence of liquid water). However, the observed subsequent replacement of 
altered magnetite-carbide-sulfide nodules by fayalite could have taken place only in the presence of 
an aqueous solution, in an asteroidal environment at temperatures below ~200°C. A major paper has 
been submitted for publication [17]. 

Silica nodules in Ca,Al-rich chondrules of the Kainsaz C03 chondrite 

In searching for Ca,Al-rich objects in the Kainsaz COS chondrite, a variety of chondrules and 
chondrule fragments with high abundance of plagioclase as well as numerous CAIs were identified in 
detailed X-ray maps of 3 thin sections. Among numerous An-rich chondrules with crystalline, silica- 
bearing mesostases, two chondrules (Fig. 1) contain unusual Si0 2 -rich nodules texturally reminiscent 
of immiscible silicate melts. 



Fig. 1. SiOrbearing chondrules in the Kainsaz C03 chondrite. White - opaque minerals, light gray - hedenbergite, 
medium gray - augite, dark gray - enstatite (En) and mesostasis, dark - a silica phase (Sil), black - voids. 


Only one of these chondrules (left panel in Fig. 1) has been studied in detail as yet. It consists 
of large subhedral grains of enstatite (20.5 vol. %) and Al,Cr-bearing augite (28.6 %) enclosed in 
plagioclase-rich mesostasis (25.8 %). Numerous nodules (24.4 %) of almost pure Si0 2 and secondary 
hedenbergite are enclosed in all phases but enstatite. Many hedenbergite nodules contain tiny relic 
grains of Si0 2 , unambiguous evidence for the substitution of hedenbergite for a silica mineral, 
probably in an asteroidal environment. An unusually Si0 2 -rich bulk chemical composition of the 
chondrule (67 wt. % Si0 2 , 5.3 % A1 2 0 3 , 1.6 % FeO, 16.5 % MgO, 6.5 % CaO, 0.9 % Na 2 0) means 
that initial crystallization of enstatite further increased the Si0 2 content of the residual melt to as 
much as -71 wt. % which is enough for a silica mineral to appear on the liquidus, but certainly not 
high enough for liquid immiscibility to occur in the CM AS system (Berman, 1983). It appears that 
enstatite and a silica phase crystallized before augite, which is consistent with the crystallization 
sequence inferred from petrographic observations, but this cannot explain the nodular appearance of 
the silica phase. High concentrations of A1 2 0 3 , CaO, Si0 2 and Na 2 0 in the chondrule imply that the 
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chondrule precursors formed in a highly fractionated nebular system, depleted in Mg and Fe. Further 
study of these unusual chondrules is in progress. 


Early solar system chronology 

Recently we began a long-term collaborative program with the Department of Earth and 
Planetary Sciences at Harvard University (Prof. S. B. Jacobsen, Dr. Q. Z. Yin) to study early solar 
system chronology. Because the distribution of 26 A1 and 53 Mn in the early solar system may have 
been inhomogeneous (e.g., Lugmair and Shukolyukov, 1998; Wood, 1998) we intend to study other 
extinct nuclide systems such as 92 Nb - 92 Zr, 182 Hf- 182 W, 97 Tc - 97 Mo, I07 Pd - 107 Ag. 

Before using the 92 Nb - 92 Zr system for precise dating of early solar processes, the initial 
abundances of both isotopes need to be measured in samples of known age. One of the best minerals 
to use for measuring Zr isotopes is meteoritic zircon. Last year we searched for zircon grains by 
making fine-scale X-ray maps of 7 polished sections of the Chaunskij mesosiderite and several thin 
sections of the Bereba, Cachari, and Pasamonte eucrites. Given the rarity and very small grain size of 
meteoritic zircon, the search for them is extremely time consuming. Eucrites contain only tiny (1-2 
|iim) zircon grains, but in Chaunskij a large (—15 pm) euhedral zircon grain was found, in addition to 
another grain that we reported earlier [18]. Laser ablation ICP-MS measurements of Zr and Nb 
isotopes in one zircon grain from Chaunskij and several rutile grains from the Zagora and Toluca 
meteorites [19] revealed that the Chaunskij zircon is 64-97 Ma older than the rutile grains. The other 
zircon grain from Chaunskij has been saved for precise U-Pb study by the SHRIMP ion probe at 
Stanford University, which is in progress. 

Wood reviewed isotopic evidence for the need to accrete planetesimals very early, while 
enough 26 A1 remained alive so that its decay could provide the heat needed to metamorphose the OC 
parent bodies and, especially, to partially melt the eucrite parent body [20]. This short time scale of 
accretion is consistent with the early time of formation of chondrules and CAIs Wood (1996) has 
argued for, during which the large amounts of energy dissipated by disk accretion would be available 
for the thermal processing of those objects. Noting the dependence of meteoriticists on outdated 
estimates of pressure and temperature profiles in the solar nebula, such as that of Lewis (1972) (28 
years old!) Wood [21] also reviewed modem astrophysical estimates of these parameters, noting that 
they changed very rapidly with early nebular evolutions and did not become quasi-stable until the 
quiescent disk phase, by which time P and T in the vicinity of the present asteroid belt were far too 
low to be consistent with chondrule and CAI formation. 
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